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Abstract 

Jet quenching as observed via high pr hadron production is studied 
in the HIJING1.35 model. We test the model on recent WA98 data on 
160 AGeV Pb + Pb ^ tv° up to 4 GeV/c. The results depend sensitively 
on the model of the Cronin effect. At (RH1C) collider energies (-y/s > 
200 AGeV), on the other hand, semi-hard hadron production becomes 
insensitive to the above model and is thus expected to be a cleaner probe 
of jet quenching signature associated with non-Abelian radiative energy 
loss. 

WA98 Pb + Pb -> 7T° data in the 2 GeV < p± < 4 GeV range were 
analyzed recently by Wang || via a parton model. He showed that the parton 
model H is capable of reproducing the observed ir° invariant inclusive cross 
sections simultaneously in p + p, S + S [|, and central Pb + Pb § in the 
CERN/SPS energy range, y/s ~ 20 AGeV. That analysis is interesting because 
it seems to imply the absence of quenching of high transverse momenta hadrons 
that should arise if partons lose energy in dense matter || |(|. The S+S and Pb+ 
Pb data clearly reveal, on the other hand, the expected Cronin enhancement JtJ 
of moderate p± pions, as seen first in p + A. The pp data confirm the need to 
supplement the LLA with non-perturbative intrinsic k±_ (^j, ^) • We address this 
problem here, following Ref. [^0), using the HIJING event generator [[jX 12 1. 
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Figure 1: WA98 results 1 on central (10%) Pb + Pb -> tt at 158 ,4GeV. Data 
(solid dots) are compared FRITIOF7.02 (open squares) 14 , VENUS4.12 (open 
circles) 13 , and parton model 2 calculations. Inset shows the ratios of the results 
of the Monte Carlo codes to data. Figure is from WA98 Coll. , Ref. 1 . 

Figure displays the prediction of VENUS ^ and FRITIOF @ models 
also for the Pb+Pb collision. It can be clearly seen that at high pt VENUS 
overpredict the experimental data of WA98 by a factor of 2.5 and FRITIOF 
underpredict it by order of magnitude. A similar discrepancy was reported for 
S+Au collisions Q. The solid curve in Figure |l| is the parton model result |2). 

We note that the spectra can also be fit well with a simple Boosted Fireball 
model j^] (see Figure ^| for S+ Au) . However such fits are far less informative 
than the parton model model one since the observed A 4 / 3 scaling of the high px 
tails is predicted and in accord with the parton model, while that dependence is 
not predicted in the fireball model. The WA98 data JlJ can also be reproduced 
by more detailed scaling hydrodynamic calculations with suitable initial and 
final freeze-out criteria [|l6|. Again the real power of pQCD is its prediction 
without parameters of the A dependence. 
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In central heavy ion reactions at the SPS, conservative estimates of the 
initial energy density suggest that e(l fm/c) » 1 ~ 5 GeV/fm 3 , are reached. 
The gluon radiative energy loss of partons in such dense matter is expected 
to exceed dE/dx ~ 1 GeV/fm Q. Recent theoretical analysis |17 predicts in 
fact a much larger non-linear energy loss AE ~ (Ax) 2 GeV/fnr if a parton 
traverses a quark-gluon plasma of thickness Ax. On the other hand, the WA98 
data seem to rule out dE/dx > 0.1 GeV/fm, as Wang showed in Ref. In 
this work, we consider this problem using the nuclear collision event generator, 
HIJING 00. 

Recall that in the LLA, pQCD predicts that the single inclusive hadron cross 
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section is given by Refs. [|[ |) 

da AB^hX 



dycPpr 



= K^2 j d 2 K a d 2 K b g(K a )g(K b ) 

abed 

dx a dx b f a/A (x a , Q 2 )f b / B {xb, Q 2 ) 



da. ^ D° h/C M 2 ) 

-f{ab^cd) . (1) 

at ttz c 

This formula convolutes the elementary pQCD parton-parton cross sections, 
da(ab — > cd), with non-perturbative lowest order fits of the parton structure 
functions, f a /A, and the parton fragmentation functions, D h / C1 to ep and e + e~ 
data. Here K a , K b are the intrinsic transverse momenta of the colliding partons. 
The model includes a K w 2 factor to simulate next-to leading order correc- 
tions p8| at a hard scale Q ~ pt/z c - The scale dependence of the structure 
and fragmentation functions account for the multiple soft collinear radiative ef- 
fects. However, below energies y/s < 100 GeV, it is well known M that LLA 
significantly underpredicts the p± < 10 GeV cross section, and additional non- 
perturbative effects must be introduced to bring LLA into agreement with data. 
Unfortunately, as emphasized in Ref. [Bj, the results are then quite model de- 
pendent below collider (RHIC) energies. As we show below, the good news is 
that this pr broadening dependence is reduced significantly at collider (RHIC) 
energies. 

In spite of the inherent ambiguity of the parton model analysis at SPS 
energies, a successful phenomenological approaches to this problem has been 
developed via the introduction || of intrinsic transverse momenta of the col- 
liding partons as in (|l]). Originally, a Gaussian form for that distribution with 
~ 1 GeV 2 was proposed in Ref. Q. However, in order to reproduce pp 
data more accurately and to take into account the Cronin effect, the Gaus- 
sian ansatz was generalized in Ref. [0 to include Q 2 and A dependence. With 
g(K) — ► g a {K,Q 2 ,A), excellent fits |g] to the WA98 data could be obtained 
assuming a factorized Gaussian distribution with 

<- 2 (b)U = (l + 0.2 Q 2 a s (Q 2 ) + °- 23fT f^ ln2Q ) (GeV/c) 2 (2) 

where Q 2 is measured here in (GeV/c) 2 . In (||) a pp tA(b) is the average number of 
inelastic scatterings a nucleon suffers traversing nucleus A at impact parameter 
b. The nuclear thickness function, t^ip), is normalized as usual to J d 2 btA(b) = 
A. Eq.(||) is the main source of the model dependence in Ref. 

The HIJING1.35 Monte Carlo model fllj, |l| incorporates pQCD jet pro- 
duction together with initial and final state radiation according to the PYTHIA 
algorithm Jl9(] . In addition, it incorporates a variant of the soft string phe- 
nomenology similar to the Lund/FRITIOF Q and DPM |HJ models to sim- 
ulate beam jet fragmentation and jet hadronization physics. Low transverse 
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momenta inelastic processes are of course highly model dependent, and the pa- 
rameters must be fit to pp and AB data O, [T^]. It is of interest to apply 
HIJING to the present study because it incorporates in addition to the above 
soft and hard dynamics, a model of soft (Cronin) multiple initial state collision 
effects as well as a simple jet quenching scheme. With these features, we are able 
to study how competing aspects of the reaction mechanism influence hadronic 
observables and explore the magnitude of theoretical uncertainties. 

In the HIJING model, excited baryon string are assumed to pick up random 
transfer momentum kicks in each inelastic scattering according to the following 
distribution 

5 (/5)(x{( K 2 +p2)( K 2 +p 2 )(1 + e ( K -p 2 )/p 3) |- 1 ; (3) 

where pi = 0.1, P2 — 1.4, p% — 0.4 GeV/c were chosen to fit low energy (p± < 1 
GeV/c) multiparticle production data in Refs. fll], Uj. A flag, IHPR2(5)(=1 or 
0) , makes it possible to compute spectra with and without this effect as shown 
in part (a) of Figure ||. The present study is the first test of this model up to 4 
GeV/c in the SPS energy range. 

Jet quenching is modeled via gluon splitting according to the number of 
mean free paths, A = 1 fm, traversed by a gluon through the cylindrical nuclear 
reaction volume. In each partonic inelastic interaction a gluon of energy AE = 
Ax dE/dx is assumed to be split off the parent jet and incorporated as a kink 
in another baryonic string The (constant) energy loss per unit length is an 
input parameter (HIPR1(14) in HIJING JT1J) and can switched on and off via 
IHPR2(4) (=1 or 0) to test the sensitivity of spectra to jet quenching as shown 
in Figure pi 

Figure^ compares the predictions of HIJING1.35 [O, |l2| without jet quench- 
ing {dE/dx — 0) for the invariant 7r° cross section in central nuclear collisions at 
SPS and RHIC energies. The cross section for central A + A collisions are com- 
puted integrating over the impact parameters up to b max chosen to reproduce 
experimental trigger cross sections. For WA98 Pb + Pb and RHIC Au + Au we 
took b max = 4.5 fm, while for the WA80 S + S we took & max = 3.4 fm. The mul- 
tiple collision eikonal geometry in HIJING is based on standard Wood-Saxon 
nuclear densities. 

The parton model fit to the WA98 data are labeled by 'Wang' from Ref. Q . 
(The normalization of both the WA98 data and Wang's latest calculations (not 
shown) have increased ~ 20 — 40% relative to Ref. ||.) We note that the 
HIJING 1.35 calculation for this interaction given by the solid jagged line also 
reproduces remarkably well the 7r° invariant cross sections without jet quench- 
ing. However, for the lighter S+S reaction, HIJING, underestimates thepj_ > 1 
GeV/c tail significantly. This error is traced to the failure of the model to repro- 
duce the pp high p± data at these energies, in contrast to its successful account 
of higher energy data [jllj . This can be seen by comparing the filled squares to 
the dot-dashed curves as explained below. Therefore, we find that the agree- 
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ment with the WA98 data is accidental and the observed A scaling of the high 
p± region at SPS energies is not reproduced. 

Figure ||a shows that the soft transverse momentum kick model is the source 
for the agreement of HIJING with the Pb + Pb data. The dot-dashed curves 
show what happens if the soft p± kicks modeled with eq.^) is turned off. The 
very strong decrease of the pion yield in the Pb + Pb case and the somewhat 
smaller but still large decrease in the S+S case shows clearly the important role 
multiple transverse kicks at these energies. In fact both S + S and Pb + Pb dot- 
dashed curves are found to coincide with the calculated pp — * tt differential cross 
section scaled by the wounded nucleon factor Wactaa/vnn, where Wa = 21, 172 
is the average number of wounded projectile nucleons and oaa = 32, 363, 636 
mb for A = 1,32,207. 

On the other hand, the data follow closely the shape of the measured pp — > 
7r + data taken from j| and scaled to AA by multiplying by the Glauber bi- 
nary collision number factor, Taa&aa- From HIJING the average number of 
binary collisions in these systems is 45, 751 resp. The fact that the WA98 data 
scale with the above Glauber factor within a factor of three, suggests that the 
additional p^ broadening due to initial state collisions is relatively small. The 
very large A dependence of the HIJING p± tail at SPS energies is due to the 
A 1 / 3 times convolution of the distribution (||). This problem is avoided in the 
parton model calculation using (0) through the separation of larger intrinsic 
momentum effects and smaller A x l dependent contributions. 

We conclude that the missing intrinsic transverse momentum component of 
HIJING precludes an accurate simultaneous reproduction of of pp, SS, PbPb 
data at SPS energies. However, unlike the parton model where no global con- 
servation laws have to be enforced, it is not clear how to incorporate intrinsic 
momenta in a global event generator like HIJING without destroying the sat- 
isfactory reproduction of low p± < 1 GeV/c data. We do not attempt to solve 
this problem here. 

Our main point is that this problem goes away fortunately at higher (RHIC) 
collider energies (y/s — 200 AGeV). As seen in Figure ||a the effect of multiple 
soft interactions is very much reduced at that energy. This is due to the well 
known effect that as the beam energy increases, the p± spectra become harder 
and additional p± smearing from initial state effects becomes relatively less 
important. It is the extreme steepness of the cross sections at SPS energies that 
amplifies so greatly the sensitivity of the moderate p± yields to this aspect of 
multiple collision dynamics. 

In Figure [|b we consider next the sensitivity of the pion yields to the sought 
after parton energy loss dynamics. The striking difference between Figure |3|a 
and Figure ||d is that in Figure ||b the SPS yield is not sensitive to the energy loss 
model in HIJING, while at RHIC energies the suppression of semi-hard hadrons 
is seen to be sensitive to jet quenching as predicted in Ref. ||. This seems 
counter intuitive at first because the increasing steepness with decreasing beam 
energy is naively expected to result in greater quenching for a fixed jet energy 
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loss. However, in this model the observed pj_ range is dominated by multiple 
soft collisions. The moderate p± < 5 GeV quarks, which fragment into the 
observed pions, are not produced in rare pQCD semi-hard interactions but are 
gently nudged several times into that p± range. Since in HIJING jet-quenching 
is restricted to only those partons that suffer a semi-hard pQCD interaction 
with a t least po — 2 GeV/c (the mini-jet scale), no quenching arises at this 
energy. 

In conclusion, we found that HIJING1.35 fits the WA98 data with or without 
jet quenching. However, this model fails to account for the weak A dependence 
of the data that scale approximately with the Glauber Taa factor, and therefore 
the fit must be viewed as accidental. The hydrodynamic model fits |l6| of WA98 
data pi unfortunately make no prediction for the A dependence. On the other 
hand, another Monte Carlo parton cascade model [0 overpredicted the ir° 
cross section at 4 GeV/c by a factor of 10. Those results and ours emphasize 
the strong model dependence of the SPS spectra due to non-perturbative aspects 
of the problem. Those aspects, while phenomenologically interesting, make it 
difficult to identify perturbative QCD phenomena and search for jet quenching. 

The main point of this work is to emphasize the good news that at collider 
energies there is much less sensitivity to the above uncertain element of the re- 
action mechanism. At > 100 AGeV the expected |(| jet quenching signature 
should be readily observable in the p± ~ 10 GeV range. Such experiments will 
commence in 1999 at RHIC and should provide important tests of the theory 
of non-Abelian multiple collisions and energy loss || |l7j ]. 



Acknowledgments 

We would like to thank T. Csorgo for the kind hospitality at Matrahaza. We also 
thank X.N. Wang, K. Eskola and D. Rischke for the extensive discussions, and T. 
Peitzmann for making the WA98 preliminary data available. This work was sup- 
ported by the Director, Office of Energy Research, Division of Nuclear Physics of 
the Office of High Energy and Nuclear Physics of the U.S. Department of Energy 
under Contract No. DE-FG02-93ER40764 and partly by US/Hungarian Science 
and Technology Joint Fund No.652/1998. and the OTKA Grant No.T025579. 

References 



[1] M.M. Aggarwal et al, WA98 Coll., nucl-ex/9806004 



[2] X.N. Wang, hep-ph/ 9804384 (1998), to be published. 



[3] J. F. Owens, Rev. Mod. Phys. 59, 465 (1987). 

[4] R. Albrecht et al, WA80 Collab., Eur. Phys. J. C5, 255 (1998) 



7 



[5] M. Gyulassy and M. Pliimer, Phys. Lett. B 243, 432 (1990); M. Gyulassy, 
M. Pliimer, M. H. Thoma and X.-N. Wang, Nucl. Phys. A 538, 37c (1992). 
M. Gyulassy and X.-N. Wang, Nucl. Phys. B 420, 583 (1994). 



[G] 
[7] 

[8] 
[9] 

[io; 
[ii 

[12 
[13 
[14 
[15 
[16 

[17 

[18 
[19 
[20 

[21 

[22 



X.-N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 1480 (1992). 

J. Cronin et al, Phys. Rev. Dll (1975) 3105; T. Ochiai, et ai, Prog. Theor. 
Phys. 75, 288 (1986); P.B. Straub et ai, E605 Collab., Phys. Rev. Lett. 68 
(1992) 452. 

R. P. Feynman, R.D. Field, and G.C. Fox, Nucl. Phys. B128, 1 (1977); 
Phys. Rev. D18, 3320 (1978). 

R. D. Field, Applications of Perturbative QCD, Frontiers in Physics Lec- 
ture, Vol. 77, Ch. 5.6 (Addison Wesley, 1989). 



M. Gyulassy and P. Levai, CU-TP-900, |hep-ph/9807247l , to be published. 

X.-N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 (1991); Comp. Phys. 
Comm. 83, 307 (1994). 

X.-N. Wang, Phys. Rep. 280, 287 (1997) 

K. Werner, Phys. Rep. 232, 87 (1993). 

B. Andersson, G. Gustafsson, and H. Pi, Z. Phys. C57, 485 (1993) 
UA. Wiedemann and U. Heinz, Phys. Rev. C56, 3265 (1997). 

A. Dumitru and D.H. Rischke, Yale preprint 1998, |nucl-th/9806003j , to be 
published. 

R. Baier, Yu.L. Dokshitzer, A.H. Mueller, S. Pcignc , D. Schiff, Nucl. Phys. 
B483 (1997) 291, B484 (1997) 265, |hep-ph/9804212| (1998) to be published. 

K. J. Eskola and X.-N. Wang, Int. J. Mod. Phys. A 10, 3071 (1995). 

T. Sjostrand, Comp. Phys. Comm. 82, 74 (1994). 

B. Andersson, et al., Nucl. Phys. B281, 289 (1987); Comp. Phys. Commun. 
43, 387 (1987) 

A. Capella, U. Sukhatme, C. I. Tan and J. Tran Thanh Van, Phys. 
Rep. 236, 225 (1994) 

K. Geiger and D. Shrivastava, Phys. Rev. C56 (1997) 2718. 



8 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

p T (GeV) p T (GeV) 

Figure 3: Invariant A + A — ► ir° cross section for central collision at SPS 
and RHIC energies are compared, a) The WA80 S + S data 4 (triangles) and 
the WA98 Pb + Pb data 1 (dots) are compared to HIJING1.35 11 with soft 
p± kicks (full lines) and without kicks (dot-dashed curves). The later scale 
with the wounded projectile number times uaa times the invariant distribution 
calculated for pp. The parton model curve is from Ref. 2 is labeled by 'Wang'. 
The filled squares show pp — > ir + data scaled by the (Glauber) number of 
binary collisions times <jaa for both SS and PbPb. b) Jet quenching, predicted 
at RHIC energies 6 , is not significant at SPS energies in the HIJING model. 
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